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Abstract 

We have investigated the escape rate of the nanospin particle with a magnetic 

field applied along the easy axis. The model studied here is described by the 

Hamiltonian H = K X S 2 + K 2 Sy + g^ b HS x , (K 1 > K 2 > 0) from which the 

escape rate is calculated within the semiclassical approximation. We have 

obtained a diagram for the orders of the phase transitions depending on the 

anisotropy constant and the external field. For K 2 /Ki > 0.85 the present 

model reveals, for the first time, the existence of the first-order transition 

within the quantum regime. 
PACS number : 75.45.+j, 75.50.Tt 
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In recent years there has been much interest in the problem of the magnetization reversal 
of a nanospin particle |T|]. It is well known that there are two possible mechanisms of the 
magnetization reversal : a classical thermal activation and a quantum tunneling 
In the study of this problem the transition rate of the system of magnetic particles is mainly 
concerned. At high temperature the transition rate is governed by the classical thermal 
activation, but when the temperature is very low the quantum tunneling dominates. As the 
temperature is lowered the phase change of the transition rate between the classical thermal 
activation and the quantum tunneling occurs. This phase transition can be either first-order 
or second-order 0. In recent several works |8|-[T0| these are shown to be possible in 
the real system such as single-domain ferromagnetic particle. These two types of phase 
transitions have been suggested by Chudnovsky fll |. However, the coexistence of the first- 



order phase transition within quantum regime and the second-order classical-to-quantum 
transition, which was also proposed by Chudnovsky, was not found yet. In this letter we 
find a real system which shows all kinds of phase transitions mentioned above. 

Consider a nanospin particle with an applied field H along the easy axis. If the spin 
particle is a uniaxial spin system with XOY easy plane anisotropy and the easy X-axis in 
the XF-plane the Hamiltonian for this system is given by || 

H = K x Sl + K 2 S 2 y - gfi b HS x (1) 

where Ki,K 2 are the anisotropy constants, fit, is the Bohr magneton, and g is the spin 
g-factor which is taken to be 2.0 here. Since we choose XF-plane as the easy plane the 
anisotropy constants satisfies Ki > K 2 > 0. 

The anisotropy energy associated with this Hamiltonian has two local energy minima ; 
the one on the + X-axis which is metastable state and the other on the —X-axis. Between 
these two energy minima there exists an energy barrier, and the spin escapes this metastable 
state either by crossing over or by tunneling through the barrier. 

In the coherent spin state representation [I2~| the effective Lagrangian corresponding to 



the Hamiltonian, Eq.(l), for small deviation from the easy plane can be written as 



c(4>,i) = ihs,-(e,4>\H\e, 



- VW ( 2 ) 



where 



m ^ = ^tFT\ \ • 2^ — oA 77 (3) 

2ivi(l — A sin — ^ cos0) 

is an effective mass, and 

V(<t>) = K 2 S 2 (sm 2 (j) + acos(f) + a) (4) 

is an effective potential for the spin system . Here we have introduced dimensionless parame- 
ters X = K 2 / Ki(< 1), a — 2H/H C (H c = being the coercive field), and added a constant 
term K 2 S 2 a in the effective potential for convenience. The effective mass and potential are 
shown in Fig.l. The barrier height of V(0) decreases as a increases and vanishes at H = H c 
for a given A. Thus, in order for the tunneling to exist a should have values between and 
2. We note that the mass depends on <p. Later we will see that the (f) dependence of the 
effective mass plays a crucial role for the occurrence of the first-order phase transition in 
quantum regime. 

At temperature T the escape rate of the spin particle can be obtained by taking ensemble 
average of the tunneling probability. Introducing Euclidean-time t = it this can be written 
as the path integral form p3Hl5| 

r(r) = J d[0(r)]e-sf W^-VbW], (5) 

where Ve{4>) = — V(<f>) is the Euclidean effective potential (see Fig.l), and = d^jdr. In 
the semiclassical approximation, neglecting the quantum fluctuation term, the escape rate 
at an energy E above the metastable minimum is given by 

T(r) ~ e-r^\ (6) 

where S(t) is the minimum effective Euclidean action which can be obtained by taking the 
smallest value of Sq and S(T). Here, So is the thermodynamic action defined by 



with E = ^r-(a + 2) 2 , and S(T) is expressed as 

r<t>f( E ) i Eh 

S(T) = 2h / d<j)\j2m(<j))[V(<j)) — E] + — , (8) 

J(j>i{E) Kb± 

where <pi{E) and 4>f(E) are the solutions of the equation V(0) = i£. For = 2aK 2 S 2 (the 
bottom of the metastable state) the Euler-Lagrange equation gives the bounce solution. 
When 2aK 2 S 2 < E < ^-(a + 2) 2 the trajectory 4>(r) in Ve(4>) shows periodic motion with 
turning points at 4>i(E) and (f)f(E). The solution corresponding to this trajectory is called 
the periodic instanton whose period is defined as 



U r4>f{E) \/2mU) 

tJE) = -^- = h dcj) / y ' . (9) 

PV ' k B T J M e) V J V {cj>)-E 



We now examine how the period r p changes as a function of energy E. Since the effective 
mass depends on cf) it influences on the variation of r p with E. To see this we look into the 
Eq.(3). For small values of A, since m(</>) varies not much it gives little effect on the behavior 
of t p . However, as A comes close to 1, the magnitude of the effective mass at turning points 
is small at first, then rises rapidly (see Fig.l). As the mass becomes larger the speed of 
a particle in the potential Ve(<P) reduces, and it takes more time to complete the periodic 
motion in Ve{4>). Thus, the period r p decreases with E at the start, but then changes to 
increase due to the rapid increase of m(<j)f). We now note that the maximum point of the 
effective mass 4> m does not coincide with the minimum point of the Euclidean potential 0o, 
which is illustrated in Fig.l. Therefore, as E approaches to the top of the potential barrier 
the motion of the particle in Ve(4>) is restricted in a region where the effective mass becomes 
small. Thus, in this region, t p decreases with E. This suggests that the whole behavior 
of t p will be the form illustrated in Fig.2. As proposed in Ref. [11] this form produces the 
first-order phase transition inside the quantum tunneling region. 

In Fig. 3 we have plotted the effective action as a function of T for a = 1 and A = 0.9. 
At T = T c the escape rate changes from the thermal activation to the quantum tunneling 



regime, and the transition is second-order. Below T c , as T is lowered, the minimum action 
increases smoothly at first, changes abruptly at T qc (the cusp at T = T qc in the figure), and 
then becomes almost constant. The quantum regime is, thus, divided into two parts : the 
thermally assisted quantum tunneling (when T qc < T < T c ) and the pure quantum tunneling 
(when T < T qc ). It can be seen from this picture that the phase transition at T = T qc is 
first-order. 

According to our numerical calculations the orders of the phase transitions are relevant 
to both a and A. In Fig. 4 we have drawn a diagram for the orders of the phase transitions in 
(A, a) plane. As remarked earlier the maximum value of a is 2 in the spin system. However, 
since we are interested in the positive effective mass a is restricted by the inequality 

A(l + ^)<1. (10) 

From the diagram we observe many interesting results. First, the classical-to-quantum 
phase transition shows both the first-order (region I) and the second-order (region II) tran- 
sitions. Note that there is only the second-order transition for A < 0.5. For materials with 
A larger than 0.5 we can see that the order changes from first to second as a increases, and 
the phase boundary increases with A up to 0.85, after which it decreases. 

In the case of the phase transition within the quantum regime there is no phase transition 
for A below 0.85. However, for A > 0.85, there exists phase transition which is first-order. We 
also observe that for 0.85 < A < 0.91 the phase boundary starts from the value corresponding 
to the maximum of the phase boundary between regions I and II and increases with A. 
When A becomes larger than 0.91, however, the phase boundary decreases with A due to 
the positive mass condition, Eq.(9). This phase boundary forms a new region III in which 
both the first-order transition within the quantum regime and the second-order classical- 
to-quantum transition coexist. Finally, in the region IV, since the negative effective mass 
begins to appear the phase transition cannot be defined. 

Our speculation about these results is as following. For materials with small A the 
height of the potential barrier is small. In an ensemble of nanospin particles each will then 
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be relaxed to reverse its magnetization easily, which leads to the smooth variation of the 
escape rate with temperature T. On the other hand, if A is large the barrier height will 
become large, and hence the spin particle will be reluctant to reverse its magnetization. In 
this case it needs the energy such as latent heat to make the spin particles ready to reverse 
their magnetizations. Therefore, the escape rate experiences the first-order phase transition. 

We can also discuss the change of the order with a for a given A. For A less than 
0.5, since the barrier height is essentially small for all values of a, the nanospin particle 
can easily reverse its magnetization, which corresponds to the second-order transition. For 
0.5 < A < 0.85, the barrier height is large at small a, but it becomes smaller as a increases. 
It is thus possible for the order to change from first to second with increasing a. We now 
consider the case A > 0.85. When a is small the barrier height is so large that only the 
high temperature first-order classical-to-quantum transition is possible. However, when a is 
moderately large, the first-order transition occurs at lower temperature, i.e., quantum region, 
with the classical-to-quantum transition being changed into second-order. This leads to the 
first-order quantum transition in quantum tunneling region. As a further increases, the 
external field makes the magnetization reversal easy (small barrier height), the situation is 
same as the case of small A. 

The Fig. 5 represents the crossover temperatures as a function of a for A = 0.9. It tells 
that both T c and T qc decrease as a gets large. This is rather obvious from the fact that the 
depth of the metastable well becomes shallow with increasing a. 

In conclusion, we have investigated the phase transition of the escape rate from 
metastable states in nanospin system with a magnetic field applied along the easy axis. 
We found the coexistence of the first-order phase transition within the quantum tunneling 
region and second-order classical-to-quantum transition for large A and a, which had not 
been observed before. Furthermore, the phase diagram for the orders of the phase transitions 
in (A, a) plane is obtained. This phase diagram can be used as a guide for the experimental 
observations. 
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FIGURES 

FIG. 1. The effective potential V ((f)) (solid), Euclidean potential Ve(4>) (dashed), and the effec- 
tive mass m{<j)) (dotted), (fro is the position at which the Euclidean potential has a minimum, while 
4> m the position at which the effective mass has a maximum. (f>i(E) and (j)f(E) are the classical 
turning points at Euclidean energy — E. 

FIG. 2. The period in Euclidean potential as a function of energy at A = 0.9 and a = 1.0, 
which shows the first-order phase transition within the quantum regime. 

FIG. 3. The actions S(T) and So as a function of temperature at A = 0.9 and a = 1. T c 
corresponds to classical-to-quantum crossover temperature and T qc to the transition temperature 
between the different quantum regime. 

FIG. 4. The phase diagram for the orders of phase transition in (A, a) plane. Region I: the 
first-order classical-to-quantum transition. Region II: the second-order classical-to-quantum tran- 
sition. Region III: the second-order classical-to-quantum transition and the first-order transition 
within quantum regime coexist. Region IV: the negative effective mass area. 

FIG. 5. Crossover temperature as a function of a. Both T c and T qc decrease with increasing a. 
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